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Magnetic fields and radio polarization of barred galaxies 
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Abstract. A three-dimensional (3D) MHD model is applied to simulate the evolution of a large-scale magnetic field 
in a barred galaxy possessing a gaseous halo extending to about 2.8 kpc above the galactic plane. As the model 
input we use a time-dependent velocity field of molecular gas resulting from self-consistent 3D N-body simulations 
of a galactic disk. We assume that the gaseous halo rotates differentially co-rotating with the disk or decreasing 
its velocity in the Z direction. The dynamo process included in the model yields the amplification of the magnetic 
field as well as the formation of field structures high above the galactic disk. The simulated magnetic fields are 
used to construct the models of a high-frequency (Faraday rotation-free) polarized radio emission that accounts 
for effects of projection and limited resolution, and is thus suitable for direct comparison with observations. 
We found that the resultant magnetic field correctly reproduces the observed structures of polarization B- vectors, 
forming coherent patterns well aligned with spiral arms and with the bar. The process initializing a wave-like 
behavior of the magnetic field, which efficiently forms magnetic maxima between the spiral arms, is demonstrated. 
The inclusion of the galactic halo constitutes a step towards a realistic model of galactic magnetic fields that 
includes as many dynamical components as needed for a realistic description. 
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1. Introduction 

Barred galaxies offer a unique opportunity to study the 
interrelations between gas flows and magnetic field struc- 
ture. The bars are known to excite strong spiral pat- 
terns in the interstellar gas and to cause strongl y non - 
axisymmetric gas flows (Englmaier & Gerhard 1997 ). 
Recent observations show that the magnetic field is 
strongly modified by gas shearing motions, although the 
field perturbations do not coincide with the suspected po- 
sition of gas velocity perturbations, as expected from sim- 
ple MHD models. Some concepts suggest that the non- 
axisymmetric gas flows in the bar may constitute the main 
mechanism amplifying the galaxy-scale magnetic fields 
(Lesch & Chiba 1997 ). Studying strong gaseous spiral pat- 
terns driven by the bar perturbation also offer a unique 
opportunity to understand the relationships between spi- 
ral arms and magnetic fields that are currently lively de- 
bated. 

The question of the gas reaction to the bar potential 
has a long history (Moss et al. |1998|, Moss et al. |1999 



and references therein). However, little has been done to 
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s tudy the evolution of magnetic field in such objects. Moss 
( 1998 ) analyzed the magnetic field evolution in barred 
galaxies, however his models are restricted to two dimen- 
sions, while the vertical structures are essential for the 
magnetic field evolution. Regular magnetic fields are not 
only passively deformed by large-scale gas motions in the 
galactic plane but their evolution is also driven by tur- 
bulent motions via the dynamo mechanism. The latter 
process generates galaxy-scale poloidal (hence also ver- 
tical) magnetic field components. Though this genuinely 
three-dimensional magnetic field structure can still be an- 
alyzed in two dimensions by applying a solenoidality con- 
dition, we find it of great interest to apply a fully 3D 
model of the field evolution. Moreover we also get an es- 
sential vertical component of the large-scale gas motion, 
which requires a completely three-dimensional treatment. 
A three-dimensional analysis of the interaction of the dy- 
namo process wi th sp iral arms has bee n per formed by 
Rohde & Elstner (|l998|) and Rohde et al. ( |l999|) . However, 
they considered only the turbulence enhancement in star- 
forming parts of spiral arms with no compression effects 
included, which is appropriate for the weak density wave 
objects (NGC 6946, Beck & Hoernes |1996[ ) they intended 
to explain. A comprehensive study of interactions of mag- 
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netic fields with bar-induced gas flows and spiral arms 
involving the dynamo process (hence carried out in three 
dimensions) is still urgently needed. 

In our previous work (Elstner et al. (2000), v. Linden et 



Table 1. Input parameters for N-body model 



al. ( 199S| ), Elstner et al. ( 1998| )) we analyzed the magnetic 
field evolution under the influence of the bar and spiral 
arms. However, these studies assumed a thin gaseous disk 
only. A large gaseous halo may be of crucial importance 
for the problem. The dynamo process includes the gener- 
ation of large-scale poloidal fields then transformed into 
azimuthal ones by differential rotation. A realistic analysis 
of the magnetic field evolution under a combined effect of 
the gas flows and the turbulent dynamo needs to involve 
an extended gaseous halo. 

For these reasons we performed the study of magnetic 
field evolution assuming a galaxy composed of a disk and 
of a rarefied ionized halo with a scale height of 2.8 kpc. 
This expands our previous results involving a thin disk 
only. As previously, we use a realistic model of the gas fiow, 
by analyzing motions of inelastically-coUiding gas clouds 
in the self-consistent calculated gravitational potential of 
stars and gas. We adopted conditions suitable for the for- 
mation of a stellar bar. We assumed the magnetic field to 
be partially coupled to the gas via the turbulent diffusion 
process. We checked the results for all the assumptions 
by varying the diffusion coefficient, the strength of the 
dynamo action, and also by switching off the turbulent 
dynamo. 

The magnetic field in external galaxies is usually stud- 
ied by means of the radio polarization, yielding the view 
of the magnetic field integrated along the line of sight and 
convolved to a certain radio telescope beam. As our mod- 
els yield complex three-dimensional magnetic field struc- 
tures, we decided to analyze our magnetic field structures 
by simulating the polarization maps, which yields much 
more concise information, easy to compare with observa- 
tions. 

2. The model 

2.1. Numerical methods 

Our experiment consists of two basic steps: the three- 
dimensional N-body, sticky particle calculations yielding 
a realistic gas velocity field and a niagnetohydrodynam- 
ical (MHD) model of the magnetic field evolution using 
the gas flow picture from the first simulations (Elstner et 
al. ^000| , 



Linden et al. 1998, Otmianowska-Mazur et al. 



1997). The model galaxy consists of two essential com- 



ponents: a collision-less self-gravitating disk composed of 
star clusters as well as highly inelastic colliding gas clouds 
moving in the gravitational potential of stars and g as, as 
well as an analytical contribution due to the bulge and 
halo. The code involves a sticky- partic le cloud scheme as 
described by Combes & Guerin ( 1985| ). In our model the 
clouds were assumed to represent the large-scale veloc- 
ity field of the gas (see Otmianowska-Mazur et al. 1997 



Model 



Mass in 

- dark-halo Mh 9.6 

- disk Ma 7.2 

- gas Mg 1.6 

- bulge Mt 4.8 

- (Md + Mg)/Mtot 0.38 
Scale length in kpc: 

- dark halo 15 

- bulge 1.2 

- disk 6 



netic field evolution by analyzing numerically the time- 
dependent solutions of the dynamo equation: 



dB/dt = rot('U x B) + rot{a B) ~ rot(?7 rotS) 

where B is the magnetic induction, v is the large-scale 
velocity of the gas, a is the dynamo coefficient (a tensor in 
general, having a diagonal form in our case) describing the 
mean helicity of turbulent motions and rj is the coefficient 
of a turbulent diffusion. 

Calculations of the magnetic field are done with a mod- 



ified version of the ZEUS3D code (Stone & Norman, |992|) 
neglecting the action of the magnetic field on the mean gas 
flow. The effects of the magnetic diffusion as well as the 
term of the a-effect are introduced dire ctly in to the code 
(von Linden et al. [l998| , Elstner et al. |2000| ). We intro- 
duced the galactic halo of ionized gas extending to 2.8 kpc 
above the galactic plane. This value was chosen due to 
computational constraints and is four times larger than 
the disk thickness used in our earlier works (von Linden 
et al. 1998, Elstner et al 



2000) 



for the detailed model description) . We modeled the mag- 



To analyze the general properties of the magnetic field 
we decided to use the simulated maps of polarized inten- 
sity, constructed on the basis of our results. We performed 
it in the way similar to that described in our previous 
works (e.g. Otmianowska-Mazur et al. 2000) namely, by 
integrating the Stokes parameters I, Q and U along the 
line of sight, convolving them with an assumed beam and 
computing the maps of polarized intensity and the angle 
of polarization B-vectors. This process yields information 
averaged over the line of sight, hence more concise, and 
directly comparable with observations. 

2.2. Model assumptions and input parameters 

The clouds span the mass range between 5 x 10^ and 
5 X 10^ M0 distributed in 10 mass-bins with logarith- 
mic intervals. The disk of the galaxy consists of 38 000 
star clusters (about 2 • 10^ M© each) and 38 096 clouds. 
The time step of calculations is 10^ years. The model pa- 
rameters adopted in the computations are summarized in 
Table |l|. 
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Fig. 1. Magnetic field and gas density in the midplane for model Al at 0.35 Gyr (left). Velocity vectors in the corotating 
frame with the bar and magnetic energy in greyscale at the midplane (right). 



Table 2. Input parameters for MHD models 



Model 


Al 


A2 


A3 


B 


C 


aoxx = aoyy [km/s]: 


10 


1.0 


0.1 





10 


aozz [km/s]: 


-10 


-1.0 


-0.1 





-10 


rjo \cvr? Is\. 


3 X lO^'' 


3 X lO^'' 


3 X 10^^ 


3 X lO^'' 


3 X lO^"^ 


rotation of halo in Z direction: 


const 


const 


const 


const 


decreasing 



The kinematic evolution of the magnetic field is solved 
in a 3D rectangular grid of points, where the XY plane is 
the galactic plane and the Z axis is the axis of rotation. 
For all experiments we use 0.3 kpc as a grid step along 
the X and Y axes and 0.08 kpc along the Z axis. The disk 
scale height is 0.8 kpc, and halo height is 2.8 kpc above 
the galactic plane. The number of grid points in X and Y 
directions is 101, while in Z, 71 points. The adopted time 
step is 0.5 X 10^ yr. 

We performed simulations of magnetic field evolution 
in a barred gas-rich galaxy with five different sets of pa- 
rameters. The input parameters for all models are sum- 
marized in Table |^. We include the diagonal components 
of the a tensor changing sinusoidally in the disk along the 
Z-direction with the maximum value of ao presented in 
Table |. In the halo we assume axx = ctyy = oizz = above 
1 kpc. In order to saturate the dynamo process at equipar- 
tition field strength Bf,q we adopt the usual a-quenching 
with a = ao/(l + {B/B^qf) (Elstner et al. 2000). The 
diffusion coefficient 77 grows linearly in the Z-direction, 
which is the result of the hydrostatic equilibriu m assu mp- 
tion on the galactic disk (Frohlich & Schultz, 1996 ). Its 
midplane value of 3 x lO^^cm^/s results from a turbulent 
velocity of 10 km/s and a correlation time of 3 x lO^ears. 



Due to uncertainty in the correlation time and its possi- 
ble reduction by the magnetic fields we assume that the 
value of 77 is ten times smaller in experiments A2 and A3 
(without 77- quenching ) . In Al-B models we apply differ- 
ential rotation which is independent of the height above 
the disk plane up to 2.8 kpc. Such an assumption is based 
on the observational results derived by different authors 
that the Galactic halo co-rotates or nearly co-rotates with 
the di sk up to 3-4 kpc above the galactic disk (Kalberla & 
Kerp 1998 and references therein). In model C the speed 
of the differential rotation in the halo decreases exponen- 
tially with increasing distance from the disk boundary to 
half of its disk value, accor ding t o recent observation of 
NGC 5775 (TaUmann et al. ^000| ). In order to study pos- 
sible changes of the magnetic field structure caused by 
the bar and spiral arms we start our simulations from 
an already-evolved magnetic field configuration computed 
using the initial velocity field. 

The polarization maps at selected time steps of our 
models are constructed with the inclination of 30° (almost 
face-on) and 85° (almost edge-on). The model assumes 
the relativistic electrons with an energy spectral index 
7 = 2.8. Their distribution is assumed to change vertically 
as a Gaussian with the scale-height of 1 kpc. Radially, the 




Fig. 2. Polarization map of model Al at time 0.35G yr and inclination of 30° superimposed onto gaseous map (grey 
plot) integrated along the line of sight. Enlargement of the most interesting region is shown on the right. 



electron density is uniform up to i? = 12 kpc. We adopt 
70% as a maximum p olariz ation degree of synchrotron 
radiation (Pacholczyk 197C, Otmianowska-Mazur et al. 



200C) and neglect any Faraday and absorption effects. The 
maps are convolved with a beam of 10" and 15" for the 
edge-on (85°) and face-on (30°) inclination, respectively. 

3. Results 

3.1. Magnetic field structure in selected planes 

The bar turns rigidly with a corotation radius of about 
8 kpc. The gaseous spiral arms arc excited by a bar per- 
turbation but they follow the bar at a somewhat lower 
angular speed. They break into parts, some of them again 
join the opposite bar end while the others slowly disperse. 
This is accompanied by a continuous creation of new arms. 
The gas density maxima, being the regions of strongest 
cloud coagulation, can be identified with dust lanes ac- 
companying spiral arms in galaxies. Shock fronts, which 
are seen in hydrodynamical gas simulations with a bar 
gravitational potential, could be identified with the gas 
density enhancements seen in our simulations, despite the 
limited resolution of our sticky particle model. The veloc- 
ity in the corotating frame shows a very similar behavior 
to the hydrodynamical models with a strong change in 
velocity direction at the location where the shock is ex- 
pected. 

All calculated models yield qualitatively similar mag- 
netic field structures, so we decided to present figures 
mainly for the model Al. The large-scale magnetic field 
quickly responds to non-axisymmetric velocity perturba- 



tions of the gas flow in the bar and spiral arms. The mag- 
netic field is locally amplified by a local compression near 
leading sides of the bar and inner edges of the spiral arms. 
It already forms magnetic arms at early stages of evolu- 
tion. As the gaseous arms dissipate, the magnetic ones re- 
main, often between newly formed gaseous arms. In Fig. ^ 
where magnetic vectors are superimposed onto the gas 
density grey plot at the time step of 3.5 x 10^ yr, we can 
see the magnetic maxima related to newly created arms 
as well as those corresponding to arm segments present 
some time ago. The magnetic arm amplified at the earlier 
stages (about 100 Myr ago) is visible in the south-west 
part of Fig. |l| (left). The pitch angles of the magnetic vec- 
tors in the interarm space are similar to those of gaseous 
arms. The time at which the results are prsented in the 
figures was chosen so to provide the best resemblance to 
observations shown in Fig. ^. 

In Fig. |l| (right) the velocity vectors (in a co-rotating 
frame) are superimposed onto the magnetic intensity grey 
plot showing the character of gaseous flow at the same 
time interval. In region near the southern end of the bar 
we can see that the radial velocity changes its sign, causing 
the amplification of the magnetic field there. Such regions 
are distributed along gaseous non-axisymmetric features 
at different time steps and they could be identified with 
gaseous shocks. 



3.2. The polarization maps 

The map of polarization vectors for a simulated almost 
face-on (30°) galaxy (model Al, see Table H), superim- 



NGC3627 8.46 GHz PI B-vect 
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posed onto the contours of polarized intensity obtained for 
the same time step (3.5x10^ yr) as in Fig. |l], is presented in 
Fig. ^ The resultant distributions of polarized vectors and 
intensity look much smoother than the vectors presented 
in Fig. |l|, which is due to beam smoothing and integration 
of magnetic field structures along the line of sight. The 
most conspicuous features in these maps are the bright 
polarized arms extending along the gaseous spirals and 
also in the interarm regions. Fig. ^ (left) presents a three- 
armed structure of polarized ridges. Two of them coincide 
with the gas density maxima constituting newly generated 
gaseous arms. Another one runs from the SW bar end to 
the interarm space towards the east and traces an older 
gaseous structure detached and dispersed 100 Myr ago. In 
this respect the magnetic arm may trace old gaseous pat- 
terns which were generated and faded away a long time 
ago in galactic history. In experiments with a high diffu- 
sion coefficient some polarized emission also fills smoothly 
the interarm space, as is often observed in spiral galaxies 
(Beck et al. [19961) . 

The polarization B-vectors are generally aligned with 
gaseous spiral arms also in the interarm space, forming a 
coherent spiral pattern with pitch angles similar to what 
we identify as the dust lanes. However, in the model Al 
the B-vectors show locally strong pitch angle differences 
between the arms and the interarm regions. This is partic- 
ularly visible for the region around R. A. 2000 ■ ll'^59'"54*'; 
Dec2ooo : 14°58'30" where the B-vectors in the arm and 
the interarm meet at almost right angles (Fig. ^ (right)). 
The interarm vectors run parallel to the bar turning sud- 
denly towards the orientation of a spiral arm, which resem- 
bles the situation in the barred galaxy NGC 3627 as pre- 
sented in Fig. |3| (Soida et al. 2001). We note particularly 
the region around R. A. 2000 ■ 

lT^0"n4^ Dec2ooo : 13°00' 
and that around R.A.2000 : 11^20"n8''; Dec2ooo : 12°59' 
in Fig. ^ where the arm and interarm B-vectors meet at 
right angles. We note also that in Fig. || around the re- 
gion of R.A.2000 : 11^59'"55^ Dec2ooo : 14°59'30", the B- 
vectors turn suddenly from the orientation perpendicular 
to the bar to that parallel to the bar. This bears some (lim- 
it ed) re semblance to the situation in NGC 1097 (Beck et 
al 1999| ). We note however that the latter galaxy shows a 
much stronger bar than that in our model. For this reason 
the validity of such a comparison is still limited and will 
be a subject of a separate study involving an analytically 
introduced bar potential. 

The bright polarized arms can be conveniently ana- 
lyzed in the frame of azimuthal angle in the disk with 
respect to the bar orientation and ln(r), in which a log- 
arithmic spiral appears as a straight line inclined by its 
pitch angle. We analyze both the position of the arms 
relative to the bar and the orientations of polarization B- 
vectors. Fig. ^ shows the time evolution of the model B at 
four time steps: 2.3, 2.5, 2.7 and 2.9 x 10^ yr. A greyplot 
of the gas density (logarithmic scale) is included, as well. 
In our frame the gaseous structures in the inner disk (thus 
more directly connected to the bar) do not change their 
position much. However, weak gas density maxima in the 
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Fig. 3. The B-vectors of length proportional to polar- 
ized intensity of NGC 3627 at 8.46 GHz superimposed 
onto the optical image by Arp (1966). A vector length of 
10"corresponds to a polarized intensity of 50 ^Jy/b.a. The 
orientations of vectors have been corrected for Faraday 
rotation effects. The map resolution is 11" (see Soida et 
alpOOll) 



outer disk (r > 8 kpc) systematically drift to the left in 
Fig. ^, being thus delayed with respect to the bar rota- 
tion. The inner parts of the magnetic arms remain tied 
to the gaseous structures and do not change the position 
much with respect to the bar. An obvious drift of mag- 
netic structures occurs in the outer disk (Fig. ^). A kind 
of "magnetic wave" move to the left, being continuously 
generated at the bar ends. 

In order to compare the time evolution of pitch angles 
of polarization B-vectors obtained for all models, we aver- 
aged the pitch angle values between the radius r = 8 kpc 
and 12 kpc along the azimuthal direction for all time steps 
(Fig. H) . For all models the absolute value of the pitch an- 
gle is quite high - from about 18° ± 2.6° for model B 
(without the dynamo mechanism) to about 22° ± 2.5° for 
our base model (Al). For the simulations with smaller 
dynamo coefficients after 300 Myr (model A2) the mean 
pitch angle grows with time from —17° to about —22° 
with a std. dev. of ±4°. 
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Fig. 4. Contour plots of polarized intensity maps and B- 
vectors from model B superimposed onto the gas density 
maps (integrated along the line of sight) shown in the 
azimuth - log(r) frame at four time steps. 
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Fig. 5. The evolution of mean magnetic pitch angle for 
models Al, A2 and B, averaged in galactocentric rings 
between 8 and 12 kpc. 



Fig. 6. The magnetic field vectors and gas density shown 
in the plane perpendicular to the galactic plane at time 
step 490 Myr. 



Fig. y presents the structure of the poloidal magnetic 
field at the disk/halo interface for the evolution time of 
490 Myr for model Al in the plane perpendicular to the 
galactic plane. Magnetic structures in the halo are formed 
quite quickly and they persist over most of the evolution 
time. The excited mode is mainly of an SO type, which 
is a consequence of the chosen form of the dynamo ten- 
sor a. In Fig. || we can see that the magnetic lines of 
force form separate loops under and above the galactic 
plane, having opposite senses of circulation. Instead of sin- 
gle, galaxy-scale loop-like structures we obtained rather 
smaller, localized structures centered at various heights 
up to 1.5 kpc above the disk plane. The growing value of 
the magnetic diffusion in the halo does not allow the for- 
mation of stronger magnetic structures higher than 2 kpc 
above the disk plane. The calculations with the decreasing 
rotational velocity (model C, see Table 2) in the halo give 
no significant differences between models having the same 
value of the other dynamo parameters (Brandenburg et 
al. 1993). 

Fig. ^ presents the edge-on view of resultant polar- 
ization B-vectors overlaid onto the polarization intensity 
contours projected almost edge-on (inclination = 85°) 
onto the sky for the model Al at the time step 510 Myr. 
The polarization map shows the vectors which are paral- 
lel to the galactic disk. This is caused by the fact that we 
see only the parts of poloidal loops closest to the plane 
where there are enough cosmic ray electrons to gener- 
ate significant synchrotron emission. The projection and 
beam smoothing effects efficiently blur all localized verti- 
cal fields leaving only its mean plane-parallel component 
seen in polarization. This resembles the majority of galax- 
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Fig. 7. The edge-on (inclination of 85°) polarized in- 
tensity map with B-vector map from model Al at time 
510 Myr. 1' corresponds to about 2 kpc in the modeled 
galaxy 
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Fig. 8. Time evolution plots of the total magnetic en- 
ergy normalized to the initial value (in logarithmic scale) 
for models Al, A3 and B. The initial field was the sat- 
urated field calculated with the primary velocity for the 
parameters of model Al. 



ies (except NGC 4631, Golla & Hummel |1994D , in which 
a disk-parallel magne tic fie ld is observed cl ose to the disk 
pla ne (D umke et al. 1995 , Dumke et al. 
al. 



19961) 



200C, Beck et 



Fig. |8| shows the magnetic energy evolution for the 
models Al, A3 and B. We started our simulations with an 
already evolved dynamo magnetic field rather than with 
a weak seed field, therefore the energy of the model Al 
remains nearly constant. The magnetic energy of model 
B with the same diffusivity decreases, because no fur- 



ther turbulent dynamo action occurs. Surprisingly, the 
magnetic energy of model A3 with a very small value 
of a and with reduced turbulent diffusivity is increasing. 
The dynamo number has the same value as in model Al. 
Therefore the magnetic energy should remain constant, 
similar to model Al. The only explanation for the growth 
is an amplification by the velocity. At that point it is dif- 
ficult to decide if the velocity itself is acting as a dynamo. 
Because of a nearly linear growth in the magnetic field 
strength (square root of the energy) it could be a tempo- 
rary increase of magnetic energy due to the shear flow. The 
vertical diffusion time for that model is about 10 Gyr. A 
consequence would be growth of the magnetic field with- 
out a major influence of a. A quenched a and r\ would 
give a further increase of magnetic energy, due to compli- 
cated 3D velocity flows. The surprise is that the pitch an- 
gle remains large even with a weak influence of a. Further 
investigation of the magnetic field evolution in a lower dif- 
fusivity medium will be presented in a forthcoming paper. 

4. Discussion 

Our experiments clearly indicate that a strong magnetic 
component associated with spiral arms persists in spite 
of introducing the large gaseous halo and of the fact that 
spiral arm perturbations occur only in a thin disk. 

The most interesting result is the fact that during 
the whole evolution time the modeled polarized intensity 
forms wave-like maxima apparently created in the com- 
presion regions in the disk, seen as the gaseous arms. 
These gaseous spirals do not corotate with the bar but 
have a somewhat smaller angular velocity. This result is 
in agreement with earlier numerical works that proposed 
spiral arms having a pattern speed smaller than t hat of 
a bar (e.g. Sellwood fc Spar ke 198S, Lind blad e t al 
Sellwood & Wilkinson |l993| , Tagger et al. |l987|) . 



1996 



The most astonishing result is that the magnetic spi- 
ral pattern rotates slightly slower than the gaseous spi- 
rals almost independent of the diffusivity. This is a com- 
pletely new result. This process causes a drift of magnetic 
field maxima from the gaseous spirals to the interarm 
space. This bears some quantitative resemblance to e.g. 
NGC 6946 (a weakly barred galaxy. Beck & Hoernes 1996| ). 
Due to this mechanism magnetic arms could work as a 
tracer of spiral arms present at earlier time stages. For the 
first time we obtained rapid turning of the magnetic pitch 
angle between arm and interarm regions, which is also ob- 
served in NGC 3627 (Soida et al. |200l| ) and in NGC 1365 
(Shoutenkov et al. 2000, Beck et al. in prep, Lindblad 
1999). This gives the impression of two separate compo- 
nents, one in spiral arms and one between them, as sug- 



gested by some observers (Soida et al. 2001). Simulations 
with different halo-to-disk mass ratios, resulting in a spi- 
ral pattern without a bar, did not show any difference 
bet ween magnetic and optical pattern speeds (Elstner et 
al. ^000|) . 

In contrast to our previous works the resultant 
magnetic field structure is fully three-dimensional with 
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poloidal loops extending into the halo. Their existence 
does not contradict observations showing disk-parallel po- 
larization vectors close to the disk plane (Dumke et al. 



1995, Dumke et al. 2000). We note however a limited appli- 



cability of our model to particular real galaxies as we did 
not introduce concrete values of parameters to reproduce 
particular observations. We restrict ourselves to mention- 
ing phenomenological features which may be worth further 
detailed studies to explain the observations in more detail. 

5. Summary and conclusions 

The evolution of large-scale galactic magnetic fields in a 
galaxy with a gaseous halo and a bar has been demon- 
strated using 3D numerical simulations. In order to solve 
the dynamo equation we used velocities obtained from 
a self-consistent N-body, sticky particle code. We found 
that: 

1. The magnetic arms persist in the presence of a large 
halo enabling a large-scale three-dimensional field. 

2. The magnetic arms are also present in interarm re- 
gions, in agreement with observations. 

3. Our calculations confirm the fact, known from earlier 
papers, that the pattern speed of the bar is higher than 
the pattern speed of spiral arms. 

4. The magnetic field may take the form of a "magnetic 
wave" with magnetic arms drifting into the interarm 
areas due to dynamical dissipation and shear amplifi- 
cation independent of the value of the turbulent diffu- 
sion. 

5. The magnetic pitch angles for all models keep the mean 
value between —15° and —25° independent of the dif- 
fusivity. 

6. The edge-on view of the modeled magnetic field agrees 
with the majority of real observations. 
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